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RESUMEN 

El objetivo principal de este estudio es determinar el desempeno que tendra 
el aislador FPT8833/1 2-1 2/8-6 utilizado en el nuevo Centro de Investigaciones 
Cientfficas de la Universidad de Fuerzas Armadas ESPE. Analisis no lineal fue 
realizado como metodo mas exacto fin determinar respuestas dinamicas de los 
aisladores. En el analisis no lineal se utilizo el modelo de tres fases propuesto por 
McVitty y Constantinou (2015) para aisladores de triple pendulo de friccion. Un set 
de apropiados siete sismos fueron seleccionado en perfiles de suelo parecido, y con 
rangos de magnitudes a las que se esperan en Quito y el Valle de los Chillos. 

Los sismos fueron escalados o normalizados utilizando el Modelo de 
Escalamiento denominado: Promedio con Pesos Ponderados utilizando varios 
periodos y modificando solo amplitudes; considerando fallas activas en una 
distancia menor a 5 km al proyecto. 



Palabras Claves: Escalamiento de sismos. Aisladores de Triple Pendulo de 
Friccion. 



ABSTRACT 

The main objective of this study is to determine the performance of the 
FPT8833/1 2-1 2/8-6 friction pendulum device used in the of the new Reserch Center 
of Universidad de Fuerzas Armadas-Espe. The nonlinear dynamic response history 
was implemented, method that represents the most accurate analysis in order to 
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determine the dynamic responses. In this analysis the model of McVitty and 
Constantinou (2015) was implemented. A set of seven appropriate ground motions 
were used, and have been registered in soil profiles and with similar magnitude 
ranges expected in Quito and Valle de Chillos. 

The ground motions were scaled or normalized using a "weighted average" 
scaling process in which several periods are considered and only the amplitude of 
the seed of motion is modified for sites within 5km from active faults to the project. 

Keyword: Scalet of earthquake. Triple Pendulum isolators friction. 



1. INTRODUCCION 



En el diseno del nuevo Centro de Investigaciones Cientfficas de la 
Universidad de Fuerzas Armadas ESPE, se consideraron dos espectros, el de 
diseno denominado DBE, que corresponde al prescrito en la Norma Ecuatoriana de 
la Construccion NEC-15, para San Rafael (z = 0.4) en un perfil de suelo“C” (Bonifaz 
201 5) y el sismo maximo considerado MCE, que se obtiene multiplicando el espectro 
DBE por 1 .4 (Constantinou etal. 2016). En la figura 1 , se presentan estos espectros. 




Period (s) 



Figura 1 Espectros de Diseno DBE y maximo sismo considerado MCE para el 
Centro de Investigaciones Cientfficas y de Post Grado de UFA-ESPE. 
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El analisis sfsmico de los diferentes bloques estructurales se realizo 
aplicando el Metodo Espectral, empleando los espectros indicados en la figura 1. 
Aguiar et at. (201 5, 201 6) 

Ahora, interesa encontrar el comportamiento no lineal del aislador friccional 
triple pendulo tipo FPT8833/1 2-1 2/8-6, con el que se esta construyendo el Proyecto, 
ante sismos impulsivos que tienen un pulso de velocidad muy alto, los mismos que 
se indicaran posteriormente. Los espectros que se obtienen con los acelerogramas 
de los sismos seleccionados tienen que ser semejantes al espectro MCE en cierto 
rango de interes para lo cual se aplica la tecnica de escalamiento de sismos. En la 
figura 2 se indica la geometrfa del aislador y los parametros utilizados para hallar la 
respuesta no lineal. 




Figura 2 Aislador FPT8833/1 2-1 2/8-6, utilizado en la construccion del nuevo 
Centro de Investigaciones Cientfficas de UFA-ESPE. 



Tabla 1 Coeficientes de roce 



Propiedad 


A 2 = As 


Ai 


A 4 


Valor Superior UB 


0.015 


0.07 


0.08 


Valor Inferior LB 


0.008 


0.05 


0.06 


Valor Nominal 


0.01 


0.06 


0.07 



En la tabla 1, se indican los coeficientes de roce del aislador FPT8833/12- 
12/8-6. Para el analisis no lineal que se presenta en este artfculo, se trabajo con los 
valores nominales, indicados en la ultima fila de la tabla 1 . 



2. SISMOS SELECCIONADOS 

Los registros que se utilizaran en el analisis no lineal de respuesta en el 
tiempo, son los presentados en la seccion 10 del reporte “LRFD-based Analysis 
and Design Procedures for Bridge Bearings in Seismic Isolators”, MCEER-1 1 -0004, 
estos terremotos fueron rotados a lo largo de la direccion paralela y normal a la falla; 
y son los indicados en la tabla 2. La seleccion de los registros sfsmicos se debe a 
las siguientes consideraciones: 
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• Los registros son de fallas cercanas a las fallas. 

• Los registros del sismo de Gazly de 1976, provienen de una region tipo 
“backward-directivity”. Los restantes registros son de regiones “forward- 
directivity”. 

• La magnitud de los registros se encuentra entre 6.7 y 7.1 . 

• Las distancias a la fuente R rup se encuentran entre 3 y 1 2 km, 

• Los registros fueron obtenidos en suelos tipo “C” y “D” 

• Los registros de Gazli de 1976 y el de Loma Prieta de 1989, son 
considerados como registros estandar con alto contenido de frecuencia. 

• Los sismos de Northridge de 1 994 y Kobe de 1 995 son sismos impulsivos y 
son considerados por su cercania a la falla. 

• Se incluye el sismo de Turquia de 1 999 porque fue altamente destructive. 



Tabla 2 Registros de sismos utilizados 



Cod. 


Nombre del sismo 


Estacion 


M w 


R 


Suelo 

Tipo 


PGV 

(cm/s) 


EQ1 


1976 Gazli, URRS 


Karakyr 


6.69 


5.46 


C 


61.50 


EQ2 


1989 Loma Prieta 


LG PC 


6.93 


3.88 


C 


77.15 


EQ3 


1989 Loma Prieta 


Saratoga-W 
Valley Coll 


6.93 


9.31 


C 


57.09 


EQ4 


1994 Northridge 


Jensen Filter 
Plan 


6.69 


5.43 


C 


73.99 


EQ5 


1994 Northridge 


Sylmar- 
Coverter Sta. 
East 


6.69 


5.19 


C 


95.07 


EQ6 


1995 Kobe, Japon 


Takarazuka 


6.90 


3.00 


D 


75.88 


EQ7 


1999 Duzce, 

Turquia 


Bolu 


7.14 


12.41 


D 


59.68 



En la tablas 3 a 5 se presentan los valores maximos de aceleracion, 
velocidad y desplazamiento registrados y el tiempo en el que se dieron esos valores ; 
la tabla 3 corresponde a los registros de los sismos de Gazli de 1976 y de Loma 
Prieta de 1989; la tabla 4 a los registros del sismo de Northridge de 1994 y la tabla 
5 contiene los parametros de los sismos de Kobe de 1995 y de Turquia de 1997. 



Tabla 3 Propiedades de los registros EQ1 , EQ2 y EQ3 



Sismo 


EQ1 (1976 Gazli) 


EQ2 (1989 L.P) 


EQ3 (1989 L.P) 


Componente 


FN 


FP 


FN 


FP 


FN 


FP 


PGA (g) 


0.60 


0.71 


0.94 


0.54 


0.40 


0.26 


Tiempo Max 
Aceleracion 
(seg) 


7.12 


8.64 


9.57 


10.37 


5.12 


6.23 


Max. Velocidad 
(cm/seg) 


64.88 


70.96 


96.84 


72.01 


71.21 


59.94 


Tiempo Max 
Velocidad 
(seg.) 


7.79 


9.90 


8.79 


7.98 


5.00 


6.13 







83 Desempeno de aislador FPT8833/1 2-1 2/8-6 empleado en la construccion... 



Max. 

Desplazamiento 

(cm.) 


24.14 


24.82 


62.62 


30.50 


20.81 


29.75 


Tiempo Max 
Desplazamiento 
(seg.) 


8.68 


9.00 


10.59 


8.34 


6.45 


6.47 



Tabla 4 Propiedades de los registros EQ4 y EQ5 



Sismo 


EQ4 (1994 Northridge- 
Jensen Filter Plant) 


EQ5 (1994 Northridge- 
Sylmar Coverter Sta East) 


Componente 


FN 


FP 


FN 


FP 


PGA (g) 


0.52 


1.07 


0.84 


0.49 


Tiempo Max 
Aceleracion 
(seg) 


7.12 


6.70 


3.50 


6.59 


Max. Velocidad 
(cm/seg) 


67.58 


64.59 


116.21 


78.23 


Tiempo Max 
Velocidad 
(seg.) 


3.21 


3.15 


3.43 


6.19 


Max. 

Desplazamiento 

(cm.) 


42.09 


21.11 


39.32 


29.19 


Tiempo Max 
Desplazamiento 
(seg.) 


2.97 


5.39 


3.02 


5.71 



Tabla 5 Propiedades de los registros EQ6 y EQ7 



Sismo 


EQ6 (1995 Kobe- 
Takarazuka) 


EQ7 (1997 Turquia-Bolu) 


Componente 


FN 


FP 


FN 


FP 


PGA (g) 


0.65 


0.70 


0.78 


0.78 


Tiempo Max 
Aceleracion 
(seg) 


6.17 


6.00 


10.76 


10.80 


Max. Velocidad 
(cm/seg) 


72.64 


83.21 


55.01 


62.63 


Tiempo Max 
Velocidad 
(seg.) 


5.27 


5.93 


11.28 


11.10 


Max. 

Desplazamiento 

(cm.) 


20.73 


26.68 


22.68 


13.55 


Tiempo Max 
Desplazamiento 
(seg.) 


5.75 


6.12 


15.96 


11.31 
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2.1 Acelerogramas seleccionados 



En las figuras 3 a 9 se presentan los acelerogramas de los sismos 
seleccionados, a la izquierda se indica la componente normal a la falla y a la derecha 
el registro de la componente paralela a la falla, se aprecia que el tiempo de duracion 
en promedio esta alrededor de los 30 seg., y que tienen pulsos muy altos, por lo que 
causaron mucho dano en las estructuras, cuando se registraron. 
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Figura 3 Acelerogramas Normal y Perpendicular a la falla del sismo de Gazli, de 
1976 en URSS. Registro: Karakyr 
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Figura 4 Acelerogramas Normal y Perpendicular a la falla del sismo de Loma 
Prietade 1989 Registro LGPC 
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Figura 5 Acelerogramas Normal y Perpendicular a la falla del sismo de Loma 
Prietade 1989 Registro Saratoga-W. Valley Coll. 
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Figura 6 Acelerogramas Normal y Perpendicular a la falla del sismo de Northridge 
de 1994 Registro: Jensen Filter Plant 
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Figura 7 Acelerogramas Normal y Perpendicular a la falla del sismo de Northridge 
de 1994 Registro: Sylmar-Coverter Sta East. 
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Figura 8 Acelerogramas Normal y Perpendicular a la falla del sismo de Kobe de 

1995 Registro: Takarazuka 
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Figura 9 Acelerogramas Normal y Perpendicular a la falla del sismo de Turquta de 

1 999. Registro: Bolu 



2.2 Espectros de respuesta elasticos 

En la figura 10, se presentan los espectros de respuesta elastica para un 
factor de amortiguamiento del 5%, que se obtuvieron con cada uno de los 
acelerogramas indicados en apartado anterior. Con Itnea de mayor espesor se 
presenta el espectro promedio, que tiene 50% de nivel de confianza. 

Ahora bien en la figura 1 1 se presentan nuevamente los espectros indicados 
en la figura 1 pero se ha anadido el espectro promedio de los 14 registros 
correspondientes a los 7 sismos seleccionados, se aprecia una buena coincidencia 
con el espectro DBE, para el rango de pertodos entre 2 y 5 segundos. Pero para el 
sismo MCE, esta bastante bajo para el rango de pertodos indicados por lo que se 
deben escalar los registros para que las ordenadas espectrales sean compatibles 
en el rango de interes seleccionado. 
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Figura 10 Espectros de respuesta elastica de los sismos considerados. 



Spectral: MCE, DBE, AVG 
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Figura 11 Espectros de diseno DBE, maximo considerado MCE y promedio de los 

sismos seleccionado. 
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3. ESCALAMIENTO DE SISMOS 

La Norma de Chile de aislacion sismica NCh2745 (2013), recomienda 
trabajar con por lo menos tres eventos stsmicos, cuyos registros deben tener 
magnitudes, distancia a la falla, fuentes del mecanismo focal y tipo de suelo en que 
fueron registrados consistentes con el espectro de diseno seleccionado. 

Por otra parte, para cada par de registros, se debe encontrar el espectro de 
respuesta elastica para 5% de amortiguamiento y combinarlos con el criterio de la 
ratz cuadrada de la suma de los cuadrados (SRSS). Los registros deben escalarse 
de tal manera que el promedio de los espectros sea mayor que 1.17 (1.3 - 0.1*1 .3) 
veces el espectro del sismo de diseno para un factor de amortiguamiento de 0.05; 
o 1.17 veces el sismo maximo considerado. Todo esto en el rango de pertodos: 
0.5 T d a 1.25 T m . Donde T D es el pertodo de la estructura para el sismo DBE y T M es 
el pertodo para el sismo MCE. 

Algo muy importante, para complementar lo indicado sobre el escalamiento 
stsmico es lo que recomienda el ASCE-7-10, referente a estructuras que se hallan 
a menos de 5 km., del piano de ruptura de una falla activa y se refiere a que los 
factores de escalamiento deben ser obtenidos unicamente con los registros de la 
componente Normal a la falla. Para otros casos se debe trabajar con el criterio SRSS 
con las dos componentes. 

4. METODO DEL PROMEDIO CON PESOS PONDERADOS 

El procedimiento denominado: “Procedimiento con pesos ponderados” 
(weighted average) utiliza la informacion sobre la aceleracion espectral en un 
numero determinado de pertodos (se trabajara con 4 pertodos). El procedimiento se 
aplica en los siguientes pasos: 



a) El primer valor de escalamiento en amplitud Fj (J=1 hasta 7) tiene la 

finalidad de disminuir el error entre la componente normal escalada y el 
Espectro MCEr, la definicion de error que esta dado en la siguiente 
ecuacion: 




( 1 ) 



Es decir, el error es la suma ponderada de los cuadrados de la diferencia 
entre la aceleracion espectral MCEr y la aceleracion espectral de la 
componente normal escalada en cuatro pertodos. Esta ecuacion da como 
resultado la siguiente expresion directa del (Factor de escalamiento): 
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W i^ MCE R (Fi )S FN (T i ) 
p — J=L 

£w,SUT,) 

i = 1 

Los periodos y factores de peso estan dados en la tabla 6. 



Tabla 6 Pesos considerados para cuatro periodos 



Periodo T/( seg) 


1.0 


2.0 


3.0 


4.0 


Peso factor w, 


0.1 


0.3 


0.3 


0.3 



b) Las componentes normales del paso anterior no cumplen con el criterio del 
ASCE 7-10. Por lo tanto se debe construir el promedio de las siete 
componentes escaladas y comparar con el espectro MCEr en el rango de 
periodos 0.5TD a 1.25TM. Calcular los valores del periodo utilizando el 
procedimiento de ELF. El espectro promedio no debe ser menor que el 
espectro MCER en el rango de los periodos de interes. Para cumplir con 
este criterio, se determina un FI factor de escalamiento, que se aplica a 
todos los movimientos. Factor final que representa el producto de los dos 
factores FjF^. 

5. SISMOS ESCALADOS PARA ESTRUCTURAS UFA-ESPE 

El periodo promedio de vibracion de los 6 bloques estructurales de la UFA- 
ESPE, que tienen aisladores sfsmicos es T D - 3.0 s, para el sismo DEB y para el 
sismo MCE es T M = 3.4 s. De tal manera que el rango de periodos de escalamiento 
se encuentra entre 1 .5 y 4.25 s., 0.5 T D - 1.25 T M . En la tabla 7 se resume el calculo 
de los factores de escalamiento. 



Tabla 7 Calculo de factores de escalamiento mediante el metodo del Promedio de 

Pesos Ponderados. 



FACTORES DE ESCALAMIENTO 




















EQ1 


EQ2 


EQ3 


EQ4 


EQ5 


EQ6 


EQ7 


Smce(* = 1) 


0,941 


0,941 


0,941 


0,941 


0,941 


0,941 


0,941 


Sfn (* = 1) 


0,815 


0,873 


0,710 


0,512 


0,732 


0,933 


0,729 


$mce (* = 2) 


0,471 


0,471 


0,471 


0,471 


0,471 


0,471 


0,471 


Sfn (* = 2) 


0,198 


0,561 


0,314 


0,411 


0,429 


0,396 


0,257 


Smce ( t= 3) 


0,314 


0,314 


0,314 


0,314 


0,314 


0,314 


0,314 


S™ (t=3) 


0,162 


0,436 


0,176 


0,421 


0,394 


0,147 


0,147 


Smce(* = 4) 


0,235 


0,235 


0,235 


0,235 


0,235 


0,235 


0,235 


Sfn ( t= 4) 


0,191 


0,326 


0,078 


0,195 


0,186 


0,071 


0,084 


fj 


1,375 


0,869 


1,461 


1,128 


1,084 


1,144 


1,528 


Fi 


1,200 


1,200 


1,200 


1,200 


1,200 


1,200 


1,200 


F,Fj = 


1,650 


1,042 


1,754 


1,353 


1,301 


1,373 


1,833 
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En la figura 12 se presentan los espectros encontrados con los registros de 
las componentes normales multiplicados por los valores de la ultima fila de la tabla 
7. Se indica ademas el espectro promedio y se observa que los espectros que se 
hallan con los registros EQ1 (Gazly 1976) y EQ7 (Turquia 1999) se encuentran por 
debajo del espectro promedio para el rango de periodos de interes. 

Ahora, en la figura 13 se presenta el espectro promedio escalado y se 
compara con el espectro MCEr, se aprecia que en el rango de interes 1 .5-4.25 s. El 
espectro promedio escalado se aproxima bastante bien al espectro MCEr. 

Por lo tanto, para hallar la respuesta no lineal del sistema de aislamiento 
cada uno de los acelerogramas se multiplicaran por el ultimo valor de Tabla 7. A 
continuacion se indica el modelo constitutivo de los aisladores FPT. 
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Figura 12 Espectros de respuesta elasticos para 5% de amortiguamiento., 
encontrados con acelerogramas escalados y espectro promedio. 
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Figura 13 Espectro promedio escalado en el rango de 1 .5 a 4.25 seg., y espectro 

MCE. 



6. MODELO DE TRES FASES 

En la figura 2 se presento, la nomenclatura del aislador de triple pendulo de 
friccion; las dimensiones allf indicadas corresponden al FPT8833/1 2-1 2/8-6, que es 
el aislador tipo del nuevo Centro de Investigaciones de la UFA-ESPE. 

Para el aislador de triple pendulo, Fenz y Constantinou (2007, 2008); Fadi y 
Constantinou (2010) presentaron un modelo de desempeno de cinco etapas. En 
base a los 12 parametros geometricos y cuatro coeficientes de roce, que definen al 
aislador. Normalmente, se tiene: 

• R 1 = R 4 » R 2 — R 3 

• Ri = R-3 < Ai = A4 

® d 4 — d 4 y d 2 — d 2 

• h 1 = h 4 y h 2 = h 3 

Con lo que se reduce el numero de parametros geometricos a 6 y los 
coeficientes de roce a 2. Esto conduce a tener un nuevo modelo de desempeno del 
aislador FTP de 3 fases. McVitty y Constantinou (201 5). Para describir el modelo de 
tres fases, se necesita presentar las siguientes ecuaciones: 
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Ri,eff = Ri - K Para i de 1 a 4. 



( 3 ) 



d* = 



Ri 



( 4 ) 



Donde R t es el radio de curvatura; h t es la altura, presentada en la figura 2; 
R ie ff es el radio de curvatura efectivo; d* es la capacidad de desplazamiento. El 
subindice i, varia de 1 a 4. A continuacion se indican las 3 fases o regimes del 
modelo. 



• Regimen I 

El desplazamiento se da solamente en las placas 2 y 3. 

0 < u < u* 



U* — 2 (/i-L — /i 2 ) f?2,e// 



( 5 ) 



F = 



W 

2 Rl.eff 



U + W 



( 6 ) 



Donde u es el desplazamiento lateral del aislador; F es la fuerza lateral 
aplicada; W es el peso que gravita sobre el aislador. A la izquierda de la figura 14 
se aprecia al aislador interior desplazandose en las superficies 2 y 3; a la derecha 
se muestra el diagrama de histeresis correspondiente. 













\ 




















Regimen 1 



Desplazamiento 

Figura 14 Desempeno de aislador en Regimen I. 
Fuente: Me Vitty y Constantinou (2015) 



• Regimen II 

El movimiento en el aislador interior Mega a los topes y empieza a 
desplazarse en las superficies 1 y 4. Normalmente es en este regimen en que 
trabaja un aislador ante un sismo de intensidad moderada y alta. Las ecuaciones se 
indican a continuacion y en la figura 15 se presenta la curva de histeresis 
correspondiente. 
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u* <u < u** 
u** — u* + 2 d{ 

W , 

F = — (u - u*) + Mi W 

z 



( 7 ) 

( 8 ) 




Desplazamiento 



Figura 15 Desempeno de aislador en Regimen II. 
Fuente: Me Vitty y Constantinou (2015) 



• Regimen III 

Este regimen se presenta cuando el sismo es extremadamente fuerte y el 
aislador interior se topa con los topes exteriores, en estas condiciones empieza a 
deslizarse el aislador interior en las superficies 2 y 3. Las ecuaciones de calculo y 
el desempeno se indica en la figura 1 6. 

u** <u < u cap 

u cap = 2d{ + 2d* 2 (9) 

F = — — (u - u") + — — (it** - it*) + Hi_W ( 1 0) 

2 * 2 , 1 // ZRieff 




Desplazamiento 



Figura 16 Desempeno de aislador en Regimen III. 
Fuente: Me Vitty y Constantinou (2015) 
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7. ANALISIS NO LINEAL 

Para el analisis no lineal se plantea un modelo bilineal, para encontrar el 
desempeno sfsmico del aislador, el mismo que no contempla degradacion de rigidez 
en la descarga ni degradacion de resistencia; muy similar al indicado a la derecha 
de la figura 15. 

El punto de fluencia del modelo bilineal, se halla con las ecuaciones que 
definen la fuerza y desplazamiento en el regimen I; y el punto ultimo se encuentra 
con la ecuacion que determina la fuerza y desplazamiento para el regimen II. Es un 
modelo conservador ya que se pudo considerar como punto ultimo las ecuaciones 
del regimen III. 

En este apartado se describe brevemente los aspectos fundamentales del 
analisis no lineal, para el modelo constitutivo bilineal. 

A la izquierda de la figura 5 se definen las 5 ramas que considera el modelo, 
en la rama uno trabaja exclusivamente en el rango elastico con una rigidez K e que 
es la pendiente de esta recta, cuando se supera la fuerza de fluencia F Y se ingresa 
a la rama dos denominada post fluencia cuya pendiente es la rigidez K p . En la rama 
dos se va controlando la velocidad, cuando esta Mega a cero se alcanza el punto de 
control 2, que se indica a la derecha de la figura 5; ingresando a la rama de descarga 
defina por 3 con una rigidez K e . 

La rama 4 es similar a la rama 2 y la rama 5 a la rama 3, pero actuando el 
sismo en sentido contrario. 

En analisis no lineal, es muy importante la determinacion de los puntos de 
quiebre, en forma exacta y para lograr este objetivo se disminuye el incremento de 
tiempo At, como se ve mas adelante. (Aguiar 201 2, Freire y Ortiz 2014) 
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Figura 17 Identificacion de ramas y puntos de quiebre del modelo bilineal. 

En la figura 18 se presenta que para el tiempo tj sobre el sistema gravita 
una fuerza F u que esta muy proxima a la fuerza de fluencia F Y , al encontrar la 
respuesta para un nuevo incremento de tiempo At, se tiene el tiempo t i+1 y se halla 
la fuerza F l+1 que es mayor a la de fluencia F Y . En este caso, se debe retornar al 
punto tj y volver a calcular la respuesta con un menor incremento de tiempo, At*. 




Figura 18 Determinacion del punto de quiebre 1 . 

El nuevo incremento de tiempo At* se obtiene por triangulos semejantes en 
funcion de las fuerzas o en funcion de los desplazamientos. 



At* = 


Fy- F; 

— ^ At 

Fj+i - Fj 


(ii) 


At* = 


40 

< 

1 

^ , 


(12) 



Ri + 1 Ri 



La aceleracion en At* se encuentra con la siguiente ecuacion, cuya 
deduccion es muy similar a la anterior. 



U* = U 1 + — — (tij+1 - U;) (13) 

At 

Donde: Uj,u i+ 1 ,u*, son las aceleraciones del suelo en los tiempos t t , t i+1 y 
en tj + At*. Normalmente se hacen varias iteracciones para llegar al punto de 
quiebre 1 y se ingresa a la rama 2. 

En la rama 2, se cambia la rigidez del sistema a K p , se determina un nuevo 
amortiguamiento c - 2 f jm K p . Donde f es el factor de amortiguamiento y m es la 
masa del sistema y se restable el valor de At. 

Como se indico en la rama 2 se controla la velocidad, estrictamente el punto 
en que cambia de signo la velocidad, ese caso significa que se paso por la velocidad 
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igual a cero, cuando esto sucede, se retorna nuevamente al tiempo t f , y se calcula 
el nuevo At*, con la siguiente expresion, que se obtiene de triangulos semejantes 
considerando que en el punto de quiebre 2 la velocidad es cero. 

At* = - — — — At (14) 

Ri — Ri+1 

Donde: qi,q i+1 , son las velocidades en los tiempos t t ,t i+1 . De igual manera 
para este menor incremento de tiempo At* se determina la aceleracion del suelo il* 
y se procede en forma similar a la determinacion del punto de quiebre 1 . 

En resumen para la determinacion de los puntos de quiebre 2 y 4 se controla 
que la velocidad sea cero y para los puntos de quiebre 3 y 5 que la fuerza del sistema 
sea F y . 



El procedimiento de calculo presentado es valido para el modelo bilineal 
indicado en la figura 1 7. 

8. METODO p DE NEWMARK 

El marco teorico del metodo /? de Newmark, se halla descrito en Aguiar 
(201 2), Chopra (2001 ). Aqui se presenta el procedimiento para el analisis no lineal. 
Partiendo de que inicialmente se halla en la rama 1; en este caso y al inicio del 
programa para t = 0, las condiciones iniciales son nulas. Con esta introduccion el 
procedimiento de calculo es el siguiente. 

i) Se determina la rigidez equivalente k 



A 1 y 

k =k H m H c 

pA t 2 P At 



(15) 



Donde: m,c,k, son la masa, amortiguamiento y rigidez del sistema; At es el 
incremento de tiempo con que se halla la respuesta; y = 0.5; /? vale 1/6 para cuando 
se considera aceleracion lineal; 1/4 para aceleracion constante y 1/8 para 
aceleracion escalonada. 



ii) Se encuentra el incremento de cargas AQ 
1 



A Q = —m ii i+l + m 



j3At 



q t + 






2 P 



kti 



J 



-C 



f ..v 4 

q,+ 



l-Z- 



V pj 



1-^- 

2/?. 



v 



a tq, 



-k q 
(16) 



iii) Se halla el incremento de desplazamiento A q. Para un sistema de un grado 
de libertad se tiene: 



( 17 ) 
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iv) Se determina el desplazamiento y velocidad en el nuevo incremento de 
tiempo. 



Ri+i — Ri+ hq 



(18) 



9 m = 



r 

P At 



Ac h + 1 + 



i~r 

, P. 



\ . f 
q,+ 



1 -+- 

2/?. 



A t q. 



V 



(19) 



v) Como se desea generalizar el calculo a las otras ramas del modelo bilineal, 
es preferible calcular la aceleracion en el nuevo incremento de tiempo 
utilizando la ecuacion diferencial del movimiento. 



Rl + 1 — 



TTl iij+i + C Cfo+i + k 
m 



( 20 ) 



Se destaca que: 



fsi+i — k 1 



( 21 ) 



f a i + 1 — c Ri+1 



( 22 ) 



Donde fs,fa son las fuerzas asociadas a la rigidez (elasticidad) y 
amortiguamiento del sistema, que van cambiando de acuerdo a la rama en que se 
hallan en el modelo bilineal. 



vi) Se determina lafuerza/Sj +1 = k q i+1 ] y se compara con lafuerzade fluencia 
F y . Si es menor significa que se continua en la rama 1 y se continua 
trabajando con la misma rigidez y amortiguamiento. Cuando la fuerza en el 
sistema es mayor que F y se determina el punto de quiebre 1, en la forma 
indicada en el apartado anterior. 

vii) Cuando ya se ha ingresado a la rama 2, se vuelve a calcular con el 
incremento de tiempo At, se cambia la rigidez y amortiguamiento, como se 
indico en el apartado anterior y se halla la respuesta en el tiempo pero esta 
vez lo que se va controlando es la velocidad. Cuando la velocidad es cero 
pasa a la rama 3, etc. La forma de calculo del punto de quiebre 2 se indico 
en el apartado anterior. 

9. RESULTADOS 

En base a la geometria del aislador, indicada en la figura 2 y considerando 
los siguientes coeficientes de roce /^ = /t 4 = 0.07; q 2 = q 3 = 0.01. Se obtuvo los 
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parametros del modelo de tres fases de Me Vitty y Constantinou (201 5) para un peso 
total que Mega al aislador W = 268 T. 

El peso indicado es el valor promedio que Mega a los aisladores. Con los 
valores de fuerzas y desplazamiento lateral, correspondientes a las dos primeras 
fases se hallo la rigidez elastica y la de post fluencia; asi como los desplazamientos 
de fluencia y ultimo (Pero asociado a la segunda fase) con lo que se tiene definido 
todo el modelo bilineal. 



Finalmente los acelerogramas considerados se multiplied por los 
coeficientes encontrados en el escalamiento de los sismos y se hallo la respuesta 
en el tiempo, la misma que se presenta a partir de la figura 1 9 hasta la figura 25. 

Antes de comentar los resultados, se debe destacar que los sismos 
considerados fueron muy fuertes y causaron mucho dano en las estructuras. A 
manera de ejemplo basta indicar que luego del sismo de Northridge de 1994, se 
creo el Comite VISION 2000 para definir la nueva filosofia de diseno sismico, en 
ese entonces para el nuevo siglo XXI. Esto debido a las grandes perdidas humanas 
y economicas que tuvieron. 

Ahora, para ver el desempeno del aislador utilizado en la construccion del 
nuevo Centro de Investigaciones Cientfficas de la UFA-ESPE, a los acelerogramas 
de esos sismos se les ha multiplicado por los factores indicados en la tabla 7, que 
se encuentran entre 1 y 2. De tal manera que los acelerogramas resultantes 
corresponden a verdaderos mega sismos y sin embargo de ello el desplazamiento 
lateral maxi mo encontrado es de 48 cm., asociado al sismo de Loma Prieta de 1 989, 
con el registro de Saratoga. 



En las figuras 19 a 25 se indica a la izquierda el comportamiento no lineal 
del aislador y a la derecha se tiene las respuestas de desplazamiento, velocidad y 
aceleracion del aislador. 



Gazli 1976 Karakyr FN 




E TIEMPO vs DESPLAZAMIENTO 





Tiempo (s) 
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Figura 19 Comportamiento del aislador y respuestas de desplazamiento, 
velocidad y aceleracion ante el sismo de Gazli 1976. Registro Karakyr. 




Figura 20 Comportamiento del aislador y respuestas de desplazamiento, 
velocidad y aceleracion ante el sismo de Loma Prieta 1989. Registro LGPC. 



Loma Prieta 1989 Saratoga FN 
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o 0.5 



0 
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Figura 21 Comportamiento del aislador y respuestas de desplazamiento, 
velocidad y aceleracion ante el sismo de Loma Prieta 1989. Registro Saratoga. 
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Figura 22 Comportamiento del aislador y respuestas de desplazamiento, 
velocidad y aceleracion ante el sismo de Northridge 1994. Registro Jensen. 




Desplazamiento (m) 



Figura 23 Comportamiento del aislador y respuestas de desplazamiento, 
velocidad y aceleracion ante el sismo de Northridge 1994. Registro Sylmar. 
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Kobe 1995 Takarazuka FN 
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Figura 24 Comportamiento del aislador y respuestas de desplazamiento, 
velocidad y aceleracion ante el sismo de Kobe 1995. Registro Takarazuka. 



Turquia 1999 Bolu FN 
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Figura 25 Comportamiento del aislador y respuestas de desplazamiento, 
velocidad y aceleracion ante el sismo de Turquia 1999. Registro Bolu. 
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El desplazamiento del aislador FPT8833/1 2-1 2/8-6, es de 56.89 cm. 
(Constantinou et at. 2015) cantidad mayor a los 48 cm., que se esperaria con el 
sismo de Loma Prieta de 1989 (Registro de Saratoga) multiplicado por 1.754. Este 
sismo es de alta frecuencia y practicamente se ha duplicado su registro, a pesar de 
ello el aislador resiste satisfactoriamente. 



10.CONCLUSIONES 

Se ha presentado el marco teorico del modelo denominado: Promedio con 
Pesos Ponderados para el escalamiento de sismos por tratarse de un modelo 
nuevo; de igual manera y con la misma razon se ha indicado el modelo de tres fases 
desarrollado por McVitty y Constantinou (2015) para definir el comportamiento no 
lineal de un aislador de triple pendulo de friccion. 

Es importante que se conozcan los modelos indicados en el parrafo anterior 
ya que representan el estado del arte, en los dos temas anotados. Unicamente para 
complementar el marco teorico se describio con bastante detalle el analisis no lineal 
de un sistema de un grado de libertad para un modelo bilineal. 

Pero el objetivo central de este articulo era conocer el desempeno del 
aislador FPT8833/1 2-1 2/8-6, que sera colocado en el nuevo Centro de 
Investigaciones de la UFA-ESPE, ante siete sismos impulsivos que causaron mucho 
dafio a nivel mundial; los mismos que han sido amplificados por los factores de 
escalamiento encontrados en el estudio, para que sean compatibles con el espectro 
maximo considerado, con esta amplificacion se convierten en mega sismos y sin 
embargo de ello la capacidad de desplazamiento del aislador FPT8833/1 2-1 2/8-6, 
es mayor a la demanda encontrada. Por lo tanto, se comportan satisfactoriamente. 
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